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Abstract 

We consider power allocations in downlink cellular wireless systems where the basestations are equipped with 
multiple transmit antennas and the mobile users are equipped with single receive antennas. Such systems can be 
modeled as multiuser MISO systems. We assume that the multi-antenna transmitters employ some fixed beamform- 
ers to transmit data, and the objective is to optimize the power allocation for different users to satisfy certain QoS 
constraints, with imperfect transmitter-side channel state information (CSI). Specifically, for MISO interference 
channels, we consider the transmit power minimization problem and the max-min SINR problem. For MISO 
broadcast channels, we consider the MSE-constrained transmit power minimization problem. All these problems 
are formulated as probability-constrained optimization problems. We make use of the Bernstein approximation to 
conservatively transform the probabilistic constraints into deterministic ones, and consequently convert the original 
stochastic optimization problems into convex optimization problems. However, the transformed problems cannot 
be straightforwardly solved using standard solver, since one of the constraints is itself an optimization problem. 
We employ the long-step logarithmic barrier cutting plane (LLBCP) algorithm to overcome difficulty. Extensive 
simulation results are provided to demonstrate the effectiveness of the proposed method, and the performance 
advantage over some existing methods. 

Index Terms 

Multiple-input single-output (MISO), interference channel, broadcast channel, power control, probability-constrained 
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I. Introduction 

The multiuser multiple-input single-output (MISO) system can be used to model a communication 
system where there is an asymmetry between the transmitters and the receivers in terms of the number 
of antennas employed. For example, in a cellular system, typically the basestation can be equipped with 
multiple antennas, whereas the mobile users are equipped with single-antennas. Then for the downlink 
transmission, when a single cell is considered, we have a MISO broadcast channel; whereas when multiple 
cells are considered, we have a MISO interference channel. Such multiuser MISO channels have been 
extensively studied in the literature. In particular, the achievable rate regions of the MISO interference 
channels have been investigated in HI [0. In (3]| (H, the problem of maximizing the sum rate of the MISO 
broadcast channel is treated. The same problem is also considered in [5], with the additional max-min 
fairness constraint. Moreover, bargaining based game-theoretic solutions are given in [|6] for a two-user 
MISO interference channel, and in for the general K-uscr case. All these works assume that perfect 
channel state information (CSI) is available at the transmitters. 

However, in practice, due to various reasons, such as estimation/quantization errors, delayed estimation, 
and limited feedback rate, the assumption of perfect CSI at the transmitter is unrealistic. In the case of 
imperfect CSI, the naive approach of treating the imperfect CSI as if it was perfect gives rise to non-robust 
design, which leads to rapid performance degradation with the increasing error levels [8j. Motivated by 
this fact, a number of efforts have addressed robust approaches that can cope with uncertainties in the 
channel knowledge. Convex optimization techniques are important tools for obtaining computationally 
efficient (exact or approximate) solutions to the robust design problems [|9]. 

Current robust design schemes can be classified into the worst-case and stochastic approaches. In the 
worst-case analysis IfTOl flTTTl Ifl2l [fT3l Ifl4ll . the channel uncertainty is considered as deterministic and norm 
bounded. The worst-case-based optimization approaches provide robustness against CSI imperfections. 
However, the actual worst case may occur with a very low probability. Hence, the worst-case approach 
may be overly pessimistic and therefore, may lead to unnecessary performance degradation. The resulting 
optimization problem sometimes does not even have a feasible solution. Even if the problem is feasible, 
the resource utilization is inefficient as most system resources must be dedicated to provide guarantees 
for the worst-case scenarios. 
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To provide more flexibility than the worst-case designs, the outage-probability-constrained robust de- 
signs have also been recently developed (HH [HU H71 [HH H91 [[20]| USD Il22]| fl23]| where the convex 
optimization tools are employed to solve the resulting stochastic optimization (also referred to as prob- 
abilistic programming) problems. In these less conservative approaches the channel state and channel 
uncertainty are considered as random processes. Compared to the worst-case approach, the stochastic 
approach achieves better average performance while keeping the probability of the worst performance low. 
Unfortunately, probability-constrained stochastic optimization are known to be computationally intractable 
except for a few special cases. In general, such optimizations are difficult to solve as their feasible sets 
are often nonconvex. In fact, finding feasible solutions to a generic probability-constrained program is 
itself a challenging research problem in the operations research community. 

In this paper, we consider several power allocation problems in multiuser MISO channels with imperfect 
CSI. We assume that the multi-antenna transmitters employ some fixed beamformers to transmit the data 
and we focus on optimizing the transmit powers of different users. In particular, for MISO interference 
channels, we treat two closely related optimization problems: one is to minimize the total transmit power 
subject to the SINR outage constraints, and the other is to maximize the achievable SINR margin under 
the power constraint. For MISO broadcast channels, we treat the MSE-constrained power minimization 
problem. All these problems are formulated as probability- constrained stochastic optimization problems. 
The major challenge in the stochastic optimization based method is to replace the probabilistic constraints 
with deterministic ones. We make use of the Bernstein approximation technique, which is a recent advance 
in the field of probability-constrained programming [E4l . to transform the probabilistic constraints into 
deterministic constraints that are conservative. The stochastic optimization problems are then transformed 
into convex optimization problems. 

The remainder of this paper is organized as follows. In Section II, we formulate the Probability- 
constrained power optimization problems in MISO interference channels. In Section III, we propose 
solutions to the stochastic optimization problems based on the Berstein approximation technique and 
the long-step logarithmic barrier cutting plane (LLBCP) algorithm. In Section IV, we treat the MSE- 
constrained stochastic power optimization problem for MISO broadcast channels. Section IV presents the 
simulation results and finally, conclusions are drawn in Section V. 



II. Power Control for MISO Interference Channels with Imperfect CSI 
A. System Model 

We consider a MISO interference channel with K transmitters and K receivers. Each transmitter employs 
M transmit antennas and each receiver is equipped with a single receive antenna. We assume that all 
receivers treat co-channel interference as noise, i.e., they make no attempt to decode the interference. 
Assuming a narrowband channel model, the received signal at receiver k is given by 

Vk = h kk x k + h kj x j + n k (1) 

m 

where x k E C M is the transmitted signal vector by the £>th transmitter, h k j E C M is the vector of complex 
channel gains between the j'-th transmitter and the k-th receiver, n k ~ M c (0,r]l) is the complex Gaussian 
noise sample. 

We assume that each transmitter employs the beamforming technique to transmit information; that is, 
we have x k = w k s k , where w k E C M is the transmit beamformer for the link between the k-th transmitter 
and the k-th receiver, and s k E C denotes the complex data symbol intended for the k-th receiver. 

In practice, the channel state information (CSI) at the receiver or transmitter is imperfect, especially for 
the transmitter-side CSI. In this paper, we assume that the transmitter only has access to the imperfect CSI 
to form the beamforming vectors. Specifically, we have the following uncertainty model for the channel 
vectors 

h k j = h k j + 8 k j (2) 



where h k j 



T 

h 1 h M 



denotes the imperfect estimate of the actual channel vector h k j, and 8 k j 



[dlj, <5^] T denotes the channel error vector, where 

5T J l -K(o,al J ). (3) 

In order to obtain robust solutions that are less sensitive to channel uncertainties, we need to explicitly 
take into account the impefect CSI. However, the mathematical problem arising from the robust beam- 
former design is in general much more complicated than the conventional non-robust design. Thus some 
simplifications are needed to make the problem tractable. In this paper, following IfTOl , |[17l , we assume 
that the transmit beamformer is of the form w k = p k g k , where p k E M + denotes transmit power, and 
9k = [dh ■■■9k I } T e denotes the unit-norm beamforming vector, i.e., ||<7 fc ||2 = 1- Here g k depends only 



on the estimated channels and it is designed in a non-robust way. For example, in the channel-matching 
approach, we have g k = h kk . On the other hand, the design of the power allocation p k is much more 
sophisticated and it depends not only on the channel estimates, but also on the channel error statistic. In 
this paper, we assume that the beam directions {g k } are fixed and focus on the power allocation problem. 

B. Problem Formulations 

We next formulate the robust power allocation problem based on the outage probability. Our goal is to 
optimize the system performance through power allocation. The system performance is usually quantified 
by its quality of service (QoS) and the resources it uses. The most common QoS metrics are the symbol 
error rate or the achievable data rate, both of which are functions of the output SINRs. The system resources 
include the transmit power and bandwidth. We consider two probability-constrained optimization problems 
as follows. 

1) SINR-constrained power minimization: The first strategy seeks to minimize the average transmit 
power subject to the QoS constraints. Given an acceptable SINR level a k and the outage probability e k 
for the fc-th transmitter and receiver pair, we aim to minimize the average transmit power while meeting 
the SINR outage constraints of all users, i.e., 

K 

mm ^Pk 

{pfe} fc =i k=l 

s.t. Pr (r fc < a k ) < e k , k = l,...,K (4) 



p k >0, k = 1, ...,K 

where Pr(A) denotes the probability of the event A, and T k denotes the SINR at the k-th receiver, given 
by 

L'/ 



\UH I 2 

Pk\h kk g k \ 



3+k 



The design parameters Sk ensures that receiver k is served with an SINR no less than a k at least (1 — 
e k ) x 100% of the time. 

2) Max-Min SINR optimization: The second strategy is to maximize the minimum SINR among all 
receivers, subject to the SINR outage constraints, and the individual transmit power constraints. We have 



the following optimization problem. 

max min a k 

{Pk,a k }k=i 

s.t. Pr (r fc < a k ) < e k , k = l,...,K (5) 
^Pk ^Pk, k = 1, 

where p fc is the given power constraint for the &-th transmitter. 

Note that the problem in © involves the individual power constraint. Alternatively we can also consider 
the total power constraint to have the following optimization problem. 

max min a k 

{Pk>a*}fcLi 

s.t. Pr (T k < a k ) ^e k , k = 1, K 

Pk > 0, fc = (6) 

k=l 

where p to t is the maximum allowable total transmit power. 
C. Background on Bernstein Approximation 

In problems ©, ©, and ©, the probabilistic constraints make the optimization highly intractable. 
The main reason is that the convexity of the feasible set corresponding to the probabilistic constraints is 
difficult to verify. To circumvent the above hurdles, we make use of the Bernstein approximation technique 
E4l to convert the probabilistic constraints to convex constraints. Next we briefly introduce the Bernstein 
approximation. 

Suppose that F(z, Q : R n x R d ->■ E is a function of z E W 1 and £ e R d . Then the probabilistic 
constraint 

Pr{F(z,C)>0}< q (7) 
can be conservatively approximated by the following 

inf t>0 {m{4)[r l F{z, c)]} - t e } < o, (8) 

or 

inf 4>0 {tlogE{^[r 1 J P(z,C)]} - tlogf?} < 0. (9) 



where ip : R — > R be a nonnegative valued, nondecreasing, convex function satisfying ip(z) > ip(0) = 1 
for any z > 0. For example, ip{z) = exp(z). Moreover, assume that for every £ G R d the function 
F(-,0 is convex. Then tE{xjj[t~ 1 F(z, £)]} — tg is convex, and thus © is convex. Indeed, since ip(-)is 
nondecreasing and convex and F(-,£)is convex, it follows that (z,t) — > ttp(t~ 1 F(z, C)) is convex. This, 
in turn, implies convexity of the expected value function E{^[t _1 F(2, ^)]}, and hence convexity of 
tE{ip[t~ 1 F(z, £)]} - tg. Similarly, t logE{^[£ -1 .F(;z, £)]} -tlogg is convex, and thus ® is convex. 

As an important special case, suppose that £ is a random vector whose components are independent 
and nonnegative. F(z, £) is affine in £, i.e., 

d 
3=1 

and the functions fj(z),j = 0, 1, d, are well defined and convex on R n . Then the probabilistic constraint 
© can be conservatively approximated by (fTTT ). 

d 

inf i>0 {tlogE{exp[t- 1 F( 2; ,C)]}-tlogf ? } = inf t>0 {/o^ + ^t logE{exp[r 1 0/ J -(2;)]}-t log g} < 0. 

i=i 

(11) 

Furthermore, the constraint in (flT) is convex. To see this, define 

d 

M/(7)=log(E[exp{ 7 o + ^07i}]) 
= lo + Y^ lo g E [ ex P(Cj7j)]- 

3=1 



Then © can be written as 



inf t>0 {tm{t~ l f{z))-t\ogg} < 0. (13) 



The function G(j,t) = £\l/(£ _1 7) — tlogg is convex in (7, t > 0) (since ^(7) is convex) and is 
nondecreasing in 70 and every 7j -, j = 1, 2, d (since Cj > 0). Since all fj(z),j = 0, 1, d are convex, 
then G(f(z),t) is convex. Due to preservation of convexity by minimization over t > 0, (fTTT) is convex. 

Note that by using the Bernstein approximation, we can convert the intractable optimization problem 
with probabilistic constraint into an explicit convex optimization. 



III. Robust Power Optimization for MISO Interference Channels 



In this section, we apply the Bernstein approximation to obtain the convex approximations to the 
probabilistic constraints in problems ©, ©, and ©, and then solve the resulting convex problems using 
the long-step logarithmic barrier cutting plan (LLBCP) algorithm. 

A. Robust SINR-constrained Power Minimization 

The major difficulty in the robust power optimization design is to convert the probabilistic constraint 
into a deterministic one. To that end we apply the Bernstein approximation to obtain the following result. 
The proof is given in Appendix A. 

Proposition 1 The following optimization problem ([141 ) is a convex conservative approximation to the 
optimization problem in ©: 

K 

s.t. 



mm 



k=l 



inf G k (p,t k ) k = l,...,K 

tk>pk 

Pk ^ 0, k = l,...,K 



(14) 



where G k (p, t k ) is defined in (|15 



G k (p,t k ) = a k r)l + t fc ^ 

j+k 

Vk\h kk g k 



h la kPj 




to 


l - t k la kP j (y 2 kj \l T g j \ 2 



-2 - log (l - t fc 1 a kPj al j \l T g 3 



2 , 7 . ,2 -*fe lo S [ 1 + h 1 Pk°kk\ lT 9k\t -*fclog(£fe) 



, and 



(15) 



A r 2 \-tT |2l 



(16) 



5. Long-step Logarithmic Barrier Cutting Plane (LLBCP) Algorithm 

Notice that the first constraint in (fT4l is itself in terms of an optimization problem. Hence, although 
the optimization problem ([141 is convex, it cannot be straightforwardly solved using a standard convex 
optimization solver. We will employ the long-step logarithmic barrier cutting plane (LLBCP) algorithm 
to solve it. 



The detailed development of the LLBCP algorithm is found in [|25l Il26]| . Here we outline the basic 
ideas of this method. Suppose that we would like to find a solution p that is feasible for ([141 ) and satisfies 
\\p — p*\\ < qj for some optimal solution p* to (PT41) . where e > is the error tolerance parameter At 
the beginning of each iteration, the feasible set, if exists, is contained in a bounded polytope. Then, we 
generate a trial point by constructing the analytic center inside the bounded polytope, and test whether 
or not the trial point belongs to the feasible set. If this trial point is not feasible, a hyperplane through 
the trial point is introduced to cut off the violated constraint(s), so that the remaining polytope contains 
the feasible set. When the trial point is feasible but not optimal, by updating the lower bound on the 
optimal objective function value of problem (fl4l) and reducing the barrier parameter, the new optimality 
constraint(s) is generated to update the polytope. Furthermore, if the hyperplanes currently in the polytope 
are deemed "unimportant" according to some criteria, they are dropped. We can then proceed to the next 
iteration with the new polytope until the termination condition is satisfied. 

Assuming that there exist the set of feasible solutions to (fl4l) . as shown in 11251 11261 . there are three 
termination conditions in the LLBCP algorithm: 

1) Termination 1: The number of hyperplanes exceeds a certain level, so that the volume of the current 
polytope would be too small to contain a small enough ball. 

2) Termination 2: The smallest slack is smaller than a certain number, so that the polytope would be 
too narrow to contain a small enough ball. 

3) Termination 3: The duality gap is enough small, so that the algorithm may be terminated with 
optimality. 

In terms of convergence, as shown in [|25l Il26l . the LLBCP algorithm terminates with a solution p 
that is feasible for ([14]) and satisfies \\p — p*\\ < e for some optimal solution p* to (fl4l) after at most 
0(K (log 2 (l/e)) 2 ) iterations, where K is the number of variables. Note that although the LLBCP algorithm 
has the same order of complexity as the algorithm in fl27l . [1281 . in practice it is computationally much 
more efficient 



'|| ■ || denotes Euclidean norm operator. 

2 It is assumed that there exist the set of feasible solutions to dl4t . and a problem dependent constant e such that 

1) The set of optimal solutions to J14t is guaranteed to be contained in the K dimensional hypercube of half-width 1/e. 

2) The set of feasible solutions to d 14b contains a full dimensional ball of radius e. 

3) It suffices to find a solution to d!4t to within an accuracy e. 
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In Fig. Q] we give a detailed flow chart of the LLBCP algorithm; and in Algorithm 1, we give the step- 
by-step description of the algorithm. The key components of the LLBCP algorithm are then elaborated 
next. 

1) Finding r-center: (Lines 5, 10, 19, 32, 41 in Algorithm 1) 

In each iteration i, we need to generate a trial point inside the poly tope V 1 = {p E ~§l K : A l p > c 1 }. 
Here, we will generate the so-called r-center of the polytope V 1 as the trial point. First we need to define 
the so-called logarithmic barrier function 

/(P,r) = S£-X>g(«») (17) 
r z — ' 

n 

where s n = a^p — c n , and is the n-th row of A. r > is the barrier parameter. For a given value of 
t\ p l (r l ) denotes the unique minimizer of this barrier function. We refer this unique point as the r-center. 
Notice that an approximate r-center is sufficient to serve as a trial point. An approximate r-center for the 
(i + l)-th iteration can be obtained from an approximate r-center for the z-th iteration by applying 0(1) 
Newton steps G9T 

2) Dropping unimportant constraints: (Lines 14-20 in Algorithm 1) 

The j-th constraint (i.e., hyperplane) is dropped only if its slack ajp — Cj has doubled since uij was 
last reset (Line 14 in Algorithm 1) and its variational quantity zuj is small (Line 16 in Algorithm 1). 
The variational quantity toj is defined as: 

&j = 2 ' J = 1 '--i N i (18) 

s j 

where N denote the number of the rows of A. The variational quantities give an indication of the relative 
importance of the constraint ajp > Cj. 

If k indexes the lower bound, 0Jk(p) — 1 (Line 12 in Algorithm 1). Otherwise, u>k(p) = afc ^ Ck (Line 
11 in Algorithm 1), where n k is initialized as ir k = k = 1, 2K, tc 2 k+i = ^VK (Line 4 in Algorithm 
1), and is updated as n k = a k p — c k (Line 22 in Algorithm 1) if m&Xk(co k (p)) > 2(Line 14 in Algorithm 
1) and zu k (p) ^0.04 (Line 16 in Algorithm 1). 

3) Checking feasibility: Given a trial point p % G V 1 , we can verify its feasibility for problem (fT4l by 
checking if it satisfies the first constraint inf G k (p l ,t k ) ^ 0,V7c. This requires solving a minimization 

tk>Pk 

problem over t k > p k . Due to the unimodality of G k (p\t k ) in t k , we can simply take a line search 
procedure to find the minimizer t* k . 

4) Cutting off the violated constraints ): (Lines 29-31 in Algorithm 1) 

If the trial point p l e V' 1 is infeasible, then a hyperplane is generated at p % as follows: 
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Algorithm 1 : LLBCP Algorithm 
l: Initialization 

2: Set e > 0, r = \, i = 0; 
3: According to (TBI) , set 

A=[J K Ik 1k] t , c=[-\1 t k -\1 t k - \s/K] T , p = oj, and a = Ap - c; 

4: Set TTfe = i, k = 1, 2if, 112K+1 = 

5: If necessary, find approximate r-center. 
6: The Iterative step 

7: if Termination 1. N > 4093iriog 2 (l/e) 

or Termination 2: min fc (a^p — c k ) < 10 _5 e 3 /[2/^ 1 - 5 log 2 (l/e)] then 



8: STOP: the best feasible point found so far is optimal. Otherwise, no feasible point is found. 
9: else 

10: Find a new approximate r-center p = p % \ 

ii: Calculate u k (p) — ak ^~° k ,Wk. 

12: Set 0Jk(p) = 1 if k indexes the lowerbound constraint that get added in Subcase 2.2. 

13: If u)j(p) > 2, calculate w-^p) as in (fT8l) . 

14: if maxk(c0k(p)) > 2 then 
15: Case 1 : 

16: if for some j, we have caj(p) < 0.04 then 

17: Subcase 1.1: 

18: Drop the hyperplane ay, 

19: Find a new approximate r-center. 

20: else 

21: Subcase 1.2: 

22: Reset 7Tj = ajp — cj, where j be an index such that uJj(p) > 2. 

23: end if 

24: end if 

25: if maxk(oJk(p)) < 2 then 
26: Case 2 : 

27: if p is not feasible in the problem (fl4l) then 

28: Subcase 2.1: 

29: For k G K, generate hyperplane(s) as in (fl~9T ). 

30: Set vr^..^ = (-Ig^) V " (^ftr) T P- 

II * drojill II ' drop II 

31: Set N^N+\kQ. 

32: Find a new approximate r-centeiH. 

33: else 

34: Subcase 2.2: 

35: if Termination 3: 1.25Nr < e then 

36: p i is the optimal solution, and STOP. 

37: else 

38: set the lower bound / = l K p — 1.25Nr on optimal objective function of (fl~4"l) . 

39: Let l prev denote previous lower bound. If l prev < I, replace l K p > l vrev by l K p > I. 

40: Set r <- Or, where 9 G (0.5, 1). 

41: Find a new approximate r-center. 

42: end if 

43: end if 

44: end if 

45: end if 
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where K = {k : G k (p\ t* k ) > 0, k — 1, K}, and Vj r fc op is the gradient of G k (p, t* k ) with respect to p 
at p l , with the k-th component given by 




5) Updating lower bound and reducing barrier parameter: (Lines 38-40 in Algorithm 1) If the point 
p % is feasible but not optimal, the lower bound / = l^p — 1.25iW to the optimal objective function 
value of problem (fl4l) is updated (Lines 38-39 in Algorithm 1), and the value of the barrier parameter 
r is reduced (Line 40 in Algorithm 1). Notice that according to the definition of (fTTT ), for a fixed value 
of r > 0, it is desirable to minimize f{p,r), leading to a balance between the objective function and 
centrality. When we need to drive the objective function value down, we just reduce the value of the 
barrier parameter r, leading to increasing emphasis on the objective function . When r is driven to zero, 
we have the convergence to an optimal solution. 

C. Robust Max-Min SINR Optimization 

We next consider the max-min SINR optimization problem in ©. Since it is difficult to verify directly 
whether problem © is convex, we use the similar method in [fT4l to solve ©. Specifically, by introducing a 
slack variable a > 0, the epigraph form of the robust max-min SINR optimization problem with individual 
power constraints © is given by 



We demonstrate that solving S(p) can be facilitated via solving a power optimization problem defined as 



which can be solved using the similar method for solving the robust power minimization problem given 
in ©. The connection between S(p) and V(p, a) is given by the following result, and the proof is given 
in Appendix B. 

Proposition 2 V(p, a) is strictly increasing and continuous in a at any strictly feasible region and is 
related to S(p) via V(p,S(p)) =1. 

Since V(p, a) is strictly increasing and continuous in a at any strictly feasible region, there exists a 
unique a* satisfying V(p,a*) = 1. It follows from Proposition 2 that solving S(p) boils down to finding 
a* that satisfies V(p,a) = 1. Due to monotonicity and continuity of V(p,a) = 1, a* can be obtained by 
a simple bi- section search. 

Finally we note that using the same approach, we can solve the problem of robust max-min SINR 
optimization problem with total power constraints given in ©. 




(21) 




(22) 
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IV. Robust Power Optimization for MISO Broadcast Channels 

In this section, we treat a related power allocation problem for a MISO broadcast system with outage 
constraints on receiver MSE. Specifically, we consider the MSE-constrained power minimization problem 
in the downlink multiuser MISO system with Gaussian channel mismatch. We adopt the Bernstein 
approximation approach to convert the probablistic constraint into a deterministic convex constraint. Note 
that the similar problem has been considered in [flD . where the Vysochanskii-Petunin inequality (VPI) is 
employed to obtain the conservative approximations to the probablistic constraints. We will demonstrate 
in Section IV the superiority of the Bernstein approach to the VPI method. 

We consider a downlink multiuser MISO system with one base station (BS) equipped with M antennas 
and K single-antenna users. The BS transmits a symbol vector s = [si, . . . , sk] t € C K , where the symbol 
s k is intended for the A;-th user. We denote the complete downlink channel as H E C KxM . The BS is 
provided only with an estimate H of H, where H has full row rank. The CSI error matrix is given by 
A = H — H, which is assumed to contain i.i.d. complex Gaussian entries, i.e., 

Skj l ~K(0,<Tlj),k= l,..,K,j = 1,...,M. (23) 
We assume that the beamforming matrix G E £, MxK [ s se t as the Moore-Penrose pseudoinverse of the 
available imperfect CSI, i.e., G = H'. Our objective is to design the diagonal power allocation matrix 

Q 1/2 = diag( A /gI, . . . , y/qic). 

— 1/2 

The A;-th user equalizes its received signal using a one-tap equalizer with coefficient q k . Thus the 
symbol estimate at the equalizer output is given by 

s k = q k 1/2 H[k, :]GQ 1/2 s + v k , k = l,...,K (24) 

where H[k, :], denoting the k-th row of H, is the A;-th user's MISO channel; v k denotes the noise sample 
at the A>th user. We assume that Ejss^} = i", and E,{vv H } = diag(r/^, . . . , r] 2 K ), where v — [vx, . . . , v K ] T . 

We use the mean-squared error (MSE) between the transmitted symbol and the receiver equalizer output 
as the QoS metric, i.e., 

MSE fc = E{|s fc -4| 2 }, k = l,...,K (25) 

We consider the following MSE-constrained power minimization problem. The objective is to minimize 
the total transmit power, subject to constraint that the probability of MSE fc being below a target value fi k 
is no less than <\> k E (0, 1). That is, 

min tr(GQG H ) 

Qto (26) 

s.t. Pr (MSE fc < fi k ) ><f) k , k = 1, K 
Now using the Bernstein approximation we can convert (1261 into a convex optimization problem, as stated 
in the following result. The proof is given in Appendix C. 

Proposition 3 The following optimization problem (|27T ) is a convex conservative approximation to the 
optimization problem in (l26T ): 

min tr(GQG H ) 

Qto (2?) 

s.t. inf Gf SE (Q, t k ) < 0, k = 1, K 

t k >o 



14 



where Gf SE (Q,t k ) is defined in © 

Gf SE (Q,t k ) 4 (^-g^ fc )-|iogdet(l-^-A fc GQG H )-t fc log(l-0 fc ), (28) 

with A fc = diag(o-| 1 , a| M ). 

We can solve problem (1271 ) using the LLBCP algorithm discussed in Section IIII-BL 

V. Simulation Results 

In this section, we present extensive simulation results to illustrate the performance of proposed 
Bernstein approximation approach to probability-constrained power optimization in wireless networks. 
First we illustrate the performance in MISO interference channels. Then we illustrate the performance in 
MISO broadcast channels and compare it with that of the VPI-based approach given in [fTTl . 

A. SINR- constrained Power Minimization in MISO Interference Channels 

We consider a MISO interference channel shown in Fig. |2] where the distance from a transmitter to the 
corresponding receiver is 200m, and the distance between the adjacent transmitters or receivers is 400m. 
The channel from the j-th transmitter to the A;-th receiver is modeled as 

200 x 3 : " 



~~ ( 1 — ) lk,jhk,j (29) 
\dk,jj 

where d k j is the distance from the j-th transmitter to the k-th receiver; 10 log 10 l k j ~ A/"(0, 8) is a real 
Gaussian random variable accounting for the large scale log-normal shadowing; h k j ~ J\f c (0 M , I M ) is a 
circularly symmetric complex Gaussian random vector accounting for Rayleigh fast fading. 
We define 

" V X 100% (30) 

Var(6™.) 

where Var(/i^) and Var(5j™) are the standard deviations of the channel fa™, and the channel error <5™-, 
respectively. For simplicity, we assume k™- = k, Vfc, j, m, and consider cases of k E {1%, 5%, 10%, 15%} 
in the following simulations. Moreover, we assume that all receivers have the same SINR level a. 

We consider the SINR-constrained power minimization problem in MISO interference channels, given 
by ©. We solve this problem by using the Bernstein approximation and the LLBCP algorithm, as 
discussed in Section III.A-B. First we consider the impact of channel error variance on the power control 
performance. Fig. |3] shows the minimum total transmit power, P™ m (a) = J2kPk( a ) versus the required 
SINR level a for the fixed outage probability e = 5% and for the different channel uncertainty levels 
k. It is seen that when the channel uncertainty increases, it takes more power to meet the SINR outage 
requirement. For a fixed channel uncertainty level k, as the target SINR value a increases, it becomes 
exceedingly difficult to meet the outage requirement; and moreover, the transmit power increases drastically 
near some limiting SINR value. This limiting value is the one which makes the optimization problem 
infeasible. Therefore the effect of imperfect CSI is more difficult to cope with when target SINR is high. 
As the channel uncertainty increases, the maximum feasible SINR value a also decreases. 
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We next consider the impact of the outage probability requirement on the power control performance. 
Fig. |4] illustrates P™ m (a) versus the target SINR level a for the fixed channel uncertainty k = 10% 
and for different outage probability values e. It is seen that as the outage requirement becomes more 
stringent, i.e., when e becomes smaller, it takes more power to meet the SINR outage requirement, and 
the maximum feasible SINR value a becomes smaller. 

B. Max-Min SINR Optimization in MISO Interference Channels 

We now consider the max-min SINR optimization problems in MISO interference channels under either 
individual or total transmit power constraint, given by dS)-©. Again the Berinstein approximation and 
the LLBCP algorithm are employed to solve the problems, as outlined in Section III.C. In Fig. [3J we 
plot the maximum achievable SINR versus the maximum allowable total transmit power for fixed outage 
probability e = 5% and for different values of the mismatched error variance k. In Fig. [6l we plot the 
maximum achievable SINR versus the maximum allowable total transmit power for fixed k = 10% and for 
different e. It is seen that for a given maximum allowable total transmit power, the maximum achievable 
SINR decreases as the channel uncertainty k increases, or as the outage probability e decreases. 

Next we consider the case where the individual transmit powers are constrained. We assume that the 
transmitters have the same maximum allowable individual transmit power p k = p total /K, where ptotai 
denotes the total transmit power. Fig. [7] shows the maximum achievable SINR versus the total transmit 
power, for fixed e = 5% and for different values of n. Fig. [8] illustrates the maximum achievable SINRs 
versus the total transmit power, for fixed k = 10% and for different outage probability e. It is seen that 
the maximum achievable SINR under the individual transmit power constraint suffers an SINR loss as 
compared to that under the total transmit power constraint. 

C. MSE-constrained Power Minimization in MISO Broadcast Channels 

In this section, we illustrate the performance of the proposed Bernstein approximation approach to 
the MSE-constrained power minimization in a downlink MISO system, as discussed in Section IV, with 
comparison with the VPI-based approach proposed in ifTTl . 

We consider a downlink MISO system with K = 3 users and the basestation is equipped with M = 3 
transmit antennas. We set a same MSE target for all users and choose //i = /i 2 = ^3 = A* from -15dB 
to -5dB. The channel coefficients are generated as i.i.d. complex Gaussian random variables with zero 
mean and unit variance. The channel error matrix is set as A[&, :] = (JkA.[k, :], k = 1,...,K, where 
A[k,:] ~K(Pm,Im). 

First, we compare the total transmit power of the Bernstein approximation approach and that of the 
VPI-based method. In Fig.[9l the minimum total transmit power against the target MSE is plotted, for three 
channel error values erj|:=1.5 x 10~ 3 , a^=10 _3 and crj?=0.5 x 10~ 3 , at the probabilistic guarantee (pk = 0.99. 
It is seen that the proposed Bernstein approximation approach outperforms the VPI-based method in 
that it achieves lower total transmit power in all cases. Additional simulations show that the Bernstein 
approximation approach significantly outperforms the VPI-based method when the channel uncertainly is 
high (i.e., larger a\), and/or the targe MSE is low (i.e., small /^), and/or the probabilistic guarantee is 



16 



stringent (i.e., high 4> k ). For instance, we can see from Fig. [9] that for cr^=1.5 x 1CT 3 , [i =-15dB, and 
4> k =0.99, the transmit power difference between the two approaches is about 7dBW. 

In order to get more insight into the behavior of these two approaches, we verify the actual probabilistic 



guarantees by plotting the histograms of the MSE in Fig. |10(a)|fT0~(b~j} The system parameters are a 2 



1.5 x 10~ 3 , fi k = — 14dB, and <p k = 0.99. The results show that under both approaches, the probabilistic 
constraints are met. In fact, the actual probability is larger than the target probability, corroborating that 
both approaches provide conservative approximations to the original probabilistic constraints. However, 
it is seen that the histogram of the VPI-based method is much more spread than that of the Bernstein 
approximation approach. In particular, the MSE realizations under the Bernstein approximation approach 
concentrate sharply around the target MSE value /x, whereas the MSE realizations under the VPI-based 
approach spread around an MSE value that is much lower than the target value /i. Hence the VPI-based 
method provides a more conservative approximation to the original probabilistic constraint, and thus needs 
more transmit power to meet the resulting more stringent MSE constraint than the actual target constraint. 

VI. Conclusions 

We have treated the problems of robust power allocation in multiuser MISO systems with imperfect 
transmitter-side CSI. The multi-antenna transmitters are assumed to employ some fixed beamformers to 
transmit data and the transmit powers need to be optimized to satisfy certain QoS constraints, taking 
into account the uncertainty in the available CSI. Specifically, for MISO interference channels, we have 
considered the transmit power minimization problem and the max-min SINR problem, subject to the 
constraints on the SINR outage probabilities. For MISO broadcast channels, we have considered the 
transmit power minimization problem subject to the constraints on the MSE outage probabilities. Our 
key contribution is to employ the Bernstein approximation to conservatively transform the probabilistic 
constraints into deterministic ones, and consequently convert the original stochastic optimization problems 
into convex optimization problems. We have provided extensive simulation results to demonstrate the 
effectiveness of the proposed robust power optimization techniques. 



A. Proof of Proposition 1 
Proof: We have 



Appendix 



Pk\h kk g k \ 2 



3+k 

M„ I 2 



< a k (31) 



®kvt + "fc p i I h kj9j I ~Pk\ h kk g k I > (32) 



Plugging © into (l32l) . we obtain (l3~3l . 



Fk (p, i k ) = akvl + "fc z2 Pj | h kj9j + $kj9j ~ Pk h kk g k + 8 kk g h 

j^k * — 



£• H 



a k r]l + ak^PjZkj - Pkikk > (33) 
i+k 
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where £ k = [£fci>£fc2> ■■■,ikK\ T ■ Note that F k (p, £ fc ) in (I33T ) is in the form of (flOT ). 

Since 5 k j,Vj are independent, £jy,Vj are independent random variables. Furthermore, from (0, it fol- 
lows that (^-g^ + fojyS'j 1 ~' A/" c (hkjSj , 0"jy 1 | 2 ) • Then by normalizing by the variance ^ | l T gj | 2 
of the real or imaginary component, we obtain the following noncentral x 2 random variable with two 
degrees of freedom 

Kkj of \ h kj9j\ 



X2 \ . T -J 



(34) 



2 | X »j| 2 l x fjl 

Thus can write the logarithm of the moment generating function of F k (p, £ fc ) as (|35 



logE[exp (t^FfcCp,^))] = logE 



*fe lQ! fc»7fc + J^logEjexp 



(T 



exp t fc + t k 1 a k Pjtkj ~ t k 1 Pkikk 

\ &k J 



X 



I exp -( 



+ logE<^ exp 



h^kvt + Yl 

i+k 



akk \-\ T a \ 2 f- l r), 
I 1 9k\ l k Vk 



2" I Mftl 



,A7,- 



h k 3 9 3 



lo § ( x - t k la kPj(y 2 kj \^ T g j \ 



h 1 Pk 


hkkdk 


2 


1 + t fc Wfc fc 


l T 9k 


2 



— I -J - log (l + f fc VfcO-fefe I l T 9fc ' 2 



(35) 



where (1351 follows from the fact that for X ~ xl W> we h ave 

exp (t^It) 1 
E {exp (tX)} = g , with t < -. (36) 

Now using (9), we obtain the Bernstein approximation (1371) to the probabilistic constraint Pr (I^ ^ a k ) ^ 
e k in ©. 



inf tfe>Pfc G fc (p, t fc ) = tfc logE [exp (t fc ^ (p, &))] - t k \oge k < 0. 



(37) 



Substituting (|35l) into ((3~7l) . we obtain (fT5T) . Note that in order to meet the condition t < ~ in (|36l) . we 
should have 



<7 



fc J 11 r„ 12,-1 



1 flfj-l t fc OfcPj < 



Hence we have 



t k > a fc max{aLp i |l T c/ | 2 } = p k . 



j^k 



(38) 

(39) 

□ 
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B. Proof of Proposition 2 

The proof follows the similar line as that in [14]. Let us denote the set of powers obtained from 
solving S(p) by {pi} and their corresponding minimal SINR by a* k . From the definition of S(p) we have 

Pk ^ Pk, V/c and min = S(p) =5> a* k ~^ S{p) 1 \/k. As a result from the definition of V(p, a) we find 

fc 

that for the choice of {pi}, the choice of b = 1 is achievable for V(p,S(p)) and therefore V(p,S(p)) ^ 1. 

Next we show that V(p,S(p)) cannot be less than one. Let us denote the set of powers obtained by 
solving V(p,S(p)) by {p* k *}. From the definition of V(p,S(p)) we clearly have a* k * ^ S(p),\/k. If 
V(p,S(p)) < 1, i.e., if max^ = c < 1, then we define the set of powers {p k } = {p* k */c}. {p k } clearly 

k P k 

satisfy the power constraints and we have their corresponding SINRs satisfying 

\h kk g k \ 2 ^t\h kk g k \ 2 



p k 



\h kk g k \ p * k *\h kk g k \ 



Since c < 1, f k > T* k *. Therefore, we have found a set of powers which satisfy the power constraints and 
yet yield a strictly larger minimal SINR compared to what the powers {p* k } obtain. This contradicts the 
optimality of {p* k *} and therefore V(p,S(p)) = 1. The strict monotonicity and continuity of V(p,a) in 
a, at any strictly feasible region, follows from a similar line of argument. □ 

C. Proof of Proposition 3 

Proof: To be consistent with ©, we rewrite the constraint Pr (MSE fc < fx k ) > <p k as Pr (MSE fc > < 
I - fa. Using d24j) and ([25]) we have 

MSE fc = q k l A[k, :}GQG H A[k, :} H + q, 1 ^. (40) 

Thus the condition MSE& > fi k becomes 

F k (Q, A[k, :]) 4 A[k, :]GQG H A[k, :} H + V 2 k - q kf i k > 0. (41) 

Note that A[k, :] is a zero-mean complex Gaussian random vector with a diagonal covariance matrix 
A fc = diag(cx^ v ...,a 2 k M ). Thus we have (@2]) Oil. 



E{ exp (t k 1 A% :}GQG H A[k, :} H ) } = 1 
1 V }) JdetfJ-^ 



det(J - f-K k GQG H ) 



(42) 

Then the Bernstein approximation to the probabilistic constraint Pr (MSE fc > fx k ) < 1 — <p k becomes 
inf tfc>0 Gf SE (Q,t k ) < 0, where Gf SE (Q,t k ) is defined in dH 

Gf SE (Q,t k ) 4 t.logE [exp (t^F k (Q,A[k,:]))]-t k log(l -fa) 

= {vl - Wk) - | logdet (i - ^A k GQG H ^j - t k log(l - fa). (43) 



□ 
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Fig. 1. Flow chart of the LLBCP algorithm. 
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Fig. 2. The simulated MISO interference channel. 
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Fig. 4. 



Minimum total transmit power versus SINR level a, for different values of e and for k 



= 10% in a MISO interference channel. 
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Fig. 10. Histograms of the MSE under the Bernstein approximation approach and the VPI-based method. = 



-14dB, and <j> k = 0.99. 



